METHOD AND APPARATUS FOR DETERMINING 

FLOW RATES 

5 Cross-Reference to Related Applications 

This application is a continuation of U.S. Serial Number 10/030,506 filed by Brown on 
9 September 2002, which is a continuation of International Application Number 
PCT/GB99/00129 filed on 6 July 2000, which claims priority firom Great Britain 
AppMcation Number 9916022.8 filed on 9 July 1999. 

10 

This invention relates to a method and apparatus for determining flow rates. In 
particular, it is concerned with the determination of the rate of flow of fluids in a conduit, 
using techniques for acquiring a distributed temperature profile in an optical fibre over a 
period of time; using this time-dependent temperature data the mass flow rates of fluids 
15 along the conduit can be determined when appropriate constants are known. These 

constants relate to a number of parameters, of which time is of particular importance, and 
also include mesisures of distance and thermal variables such as temperature, 
conductivity and specific heat. 

20 Mass flow rate .information is a very important tool for the efficient management of oil 
wells and the like. It is of course important to have reliable production data, as soon as 
possible, not only for its own sake. If flow rate data is promptiy available, it may also be 
actively used to adjust or improve the flow rates, to diagnose immediate or potential 
problems, or to trigger alarms. Significant variations in flow can be met with an 

25 appropriate management response. 
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It has been known in principle for over 25 years that thermal data can be used to derive 
mass flow rate information, and that this information is appUcable to oil field operations 

« 

and the like. Reference is made to the paper ''Use of the Temperature Log for 
5 Determining Flow Rates in Producing Wells'\ Curtis M R and Witterholt E J, Society of 
Petroleum Engineers of AIME Paper No SPE 4637, Curtis and Witterholt describe a 
method for calculating mass flow rate of a fluid up a well bore as a function of the 
temperature profile based on algorithms developed by Ramey, published as ''Well-bore 
Heat Transmission'', Ramey H, J.PetTecL, April 1962. 
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Nevertheless, as a practical matter, in the economically and commercially important field 
of determining the mass flow rate of produced or injected fluids within a well bore, 
downhole measurements have typically been made using either spinner or venturi 
techniques in one or a plxirality of locations within the production tubing. The 
1 5 equipment or devices that have been used have been either permanently installed in the 
well bore or conveyed into a measuring location by wireline. 

These currentiy used devices do however have well known disadvantages. The spinner 
device is typically run on a wireline. Use of this technique commonly involves shutting 
20 in the well for extended periods while setting up the equipment, and then running the 
sensor and cable in the well, which presents a hazard to the integrity of the well. 
Surveys of this kind are carried out infrequently, and only provide an instantaneous 
picture of the flow characteristics of the well. 

25 In order to obtain continuous flow information, it is necessary to use downhole 
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instrumentation that is permanently installed, A particular benefit of permanent 
instrumentation is that it enables a producing well to be better controlled. Venturi 
techniques, in which the pressure drop across a known orifice is measured, enable flow 
rates to be permanently monitored, but do however have limitations. Firstly, the orifice 
5 device restricts the internal diameter of the tubing. Secondly, the device relies upon two 
independent high accuracy pressure sensors, but the output of such devices has a 
tendency to drift with time. Thirdly, the venturi device must be routinely calibrated to a 
fixed fluid mix density, to ensure continued accuracy of measurement. 

1 0 For the foregoing reasons, among others, there is and has been for a long time a 

continuing need to find and develop improved methods for downhole mass flow rate 
monitoring. 

We have developed sensing and measuring equipment based on opto-electronic systems 
15 at a surface location operatively connected to fibre optic sensors deployed downhole. 
Using such systems, it is not necessary to have any electronics downhole and the fibre 
sensors can provide temperature and pressure information, while being resistant to 
temperatures up to 250 QC and above. 

20 It has been.known for over 15 years that optical fibres can report temperature 

distributions. Seefor example GB 2122337 and EP 0213872. We have now found that 
it is possible to combine, in a usefiil and practical and advantageous manner, the 
derivation of mass flow rates firom thermal data with the acquisition of thermal data by 
means of a fibre optic sensor. 

25 
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Typically, the thermal data is acquired as follows. A laser light pulse is sent down an 
optical fibre wave guide. As the pulse of hght travels along the wave guide, the thermal 
molecular vibration at each point along the length of the wave guide causes a very weak 
5 reflected signal to travel back up the fibre towards the source. An optical coupler splits 
the reflected light away firom the fibre and takes it to a detector. The time lapse between 
the launch of the light pulse and detection allows the distance of the reflection point 
down the optical fibre to be calculated, since the speed of hght in the fibre is constant 
and is known. The amplitude of the returned light is a function of the molecular 
1 0 vibration at the reflection point, increasing with increasing temperature. As reflected 

light is detected over a time period corresponding for the time taken for the light pulse to 
travel the length of the optical fibre and back, the output of the detector is effectively a 
distributed temperature profile along the whole length of the fibre. 

1 5 The present invention addresses the deficiencies of the prior art methods of determining 
downhole flow rates and provides methods and apparatus utilising a distributed fibre 
optic sensor. We have found that a single optical sensor system can provide sufficient 
thermal information to determine the mass flow rates of produced fluids within a well 
bore, using an optical fibre placed within or adjacent to the well bore, almost 

20 instantaneously; at any time, substantially continuously if required, without interference 
with production or prejudicing the integrity of the well. 

The present invention concerns aspects of the method and apparatus described below. 
The scope of the invoition extends to all novel aspects thereof whether individually or in 
25 combination with any of the other features disclosed herein. 
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More specifically, in one aspect of the invention a method of determining mass flow 
rates of fluid in a conduit located in a heat sink differing in temperature from the fluid 
may comprise obtaining a distributed temperature profile of fluid flowing along a length 
of conduit by means of optical data obtained from a length of optical fibre in thermal 
contact therewith, obtaining a profile of the heat sink temperature extemal to the conduit, 
and deriving mass flow rates of fluids in the conduit from the said profiles and from 
measured fliermal transfer parameters. 

Correspondingly, apparatus for determining mass flow rates of fluid in a conduit located 
in a heat sink differing in temperature from the fluid may comprise a length of optical 
fibre in thermal contact with the fluid, means for'^obtaining a distributed temperature 
profile of fluid flowing along a length of conduit by means of optical data obtained from 
said length of optical fibre, and means for deriving mass flow rates of fluids in the , 
conduit from the said distributed temperature profile, from a profile of the heat sink 
temperature extemal to the conduit, and from measured thermal transfer parameters. 

In a further aspect of the invention there is provided a method of monitoring the mass 
flow rates of fluids flowing in variable quantities along a length of underground conduit, 
including monitoring the said rates during both a calibration period and an observation 
period (which may include some or all of the calibration period); which method 
comprises; 

(a) establishing distributed temperature measuring apparatus comprising an optical 
fibre extending along the said length of conduit in thermal contact with the fluid 
and/or with the conduit, together with means for passing light along the optical 



fibre in the said length, receiving light emergent therefrom and interpreting 
temperature- and location-related characteristics of said emergent light in terms 
of the temperature profile of the optical fibre at a series of locations along the 
said length of conduit; 

(b) determining the natural geothermal profile along the length of the conduit to be 
monitored (the natural profile being the profile in the absence of heating or 
cooling from the conduit); 

(c) running fluid to be monitored through the said length of conduit; 

(d) in the calibration period: 

(i) measuring the actual mass flow rates of the fluid during time intervals; 

(ii) dxiring those intervals, obtaining distributed temperature profiles of the 
length of conduit by passing light along the optical fibre and interpreting 
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temperature- and location-related characteristics of light emergent from 
the said length; 

(iii) correlating the distributed temperature profiles observed in (i) with the 
flow rates measured in (ii) whereby to obtain calibration data which 
calibrate the temperature measuring apparatus in terms of mass flow rate; 

(e) in the observation period: monitoring the distributed temperature profile of the 
length of conduit by means of the distributed temperature measuring apparatus 
and jobtsiining therefrom the flow rates of the fluid in the length of conduit using 
the calibration data obtained in the calibration period. 

Correspondingly, apparatus for monitoring the mass flow rates of fluids flowing in 

variable quantities along a length of undergroimd conduit may comprise: 

(a) distributed temperature measuring apparatus comprising an optical fibre 



extending along the said length of conduit in thermal contact with the fluid and/or 
with the conduit, together with means for passing hght along the optical fibre in 
the said length, receiving light emergent therefirom and interpreting temperature- 
and location-related characteristics of said emergent light in terms of the 
temperature profile of die optical fibre at a series of locations along the said 
length of conduit; 

(b) means for determining the natural geothermal profile along the length of the 
conduit to be monitored (the natural profile being the profile in the absence of 
heating or cooling from the conduit); 

(c) optionally, means for measuring the actual mass flow rates of the fluid during 
time intervals in a calibration period; 

(d) means for correlating distributed temperature profiles obtained from (a) with 
actual flow rates whereby to obtain calibration data which caUbrate the 
temperature measuring apparatus in terms of mass flow rate; 

(e) means for monitoring the distributed temperature profile of the length of conduit 
by means of the distributed temperature measuring apparatus and obtaining 
therefrom the flow rates of the fluid in the length of conduit using the calibration 
data obtained from (d). 

Various preferred and optional features of the invention will become apparent from the 
following description. 

In an embodiment of the invention, a fibre optic distributed temperature sensor is 
installed in a well bore inside a thermally conductive tube which is suitably clamped or 
bonded to a substantial continuous fixed structure extending over the length of well bore 



in which the optical fibre is operatively deployed. Suitably, the fixed structure may be 
the conduit for the fluid whose mass flow rate is to be measured. This conduit may be 
the oil well casing, or the production tubing, or aay other similar conduit appropriate to 
the particular downhole environment where the flow rate is required to be known. 
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The tube may be filled with thermally conductive liquid, to ensure functional thermal 
contact between the optical fibre and the conduit concerned, and preferably between the 
fibre and the fluid whose flow rate is to be determined. 

1 0 The mass flow rate of the fluid is determined by the application of predetermined 

algorithms to the distributed temperature profile that is determined for the fluid. Typical 
applications of the invention are the calculation of mass flow rate in producing oil, water 
and gas wells, in a variety of different fluid combinations, and in injecting water wells. 

1 5 Among the different possible methods of determining mass flow rate firom temperature, 
two are preferred in the practice of the present invention. 

In a first preferred method, flow rate data is derived fi:om the thermal behaviour of fluids 
flowing through massive underground formations, which act as heat sinks at their natural 

20 temperatures, that is. to say their temperatures in the absence of any flow of heating or 
cooling fluid through the conduit that runs through these formations. Generally 
speaking, in a vertical well, the temperature rises more or less linearly with depth, and it 
is normally sufficiently accurate for the purposes of the present invention to treat the 
resulting geothermal gradient as being linear. As fluid flows along the conduit it is 

25 heated or cooled by conduction. The temperature at any point depends on the thermal 
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properties of the flowing fluid, of the installed completion (the production tubing and 
associated hardware within the lined well bore of a well, including such equipment as 
down-hole safety valves, packers and circulating valves) and of the surrounding 
formation, and is dependent upon flo^y rate, pressure, volume, temperature (PVT), Joxxle- 
5 Thompson effects, and frictional losses, and can be time dependent. 

It is observed that, starting from a point of interest deep within the well, which may be a 
point at which fluid in temperature equilibrium with the surrounding formation is 

■ 

introduced into the well bore, the temperature profile of the fluid above that point (as the 
1 0 fluid flows upwards through a heat sink of progressively lower temperature) takes the 
form of a curve which approaches an asymptote to a straight line parallel to the 
geothennal profile, ie of the same geothermal gradient but displaced by a certain 
temperature. The actual shape of the asymptotic curve is determined by the thermal 
properties of the system, mass flow rates, fiiction, Joule-Thompson and PVT properties 
15 of the flowing fluid as has been described in publications such as Curtis and Witterholt, 
SPE 4637, mentioned above, to which reference should be made for fiirther details. 

Appropriate algorithms for this first preferred method are: 

20 T(z,t) = Tge+GgZ-GgA+(Tfe-Tge+GgA)e'^'^ 

a 

where 

A = QD{Cf(kh+rciUf(t))/2DrciUkh 

25 
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and 



f(t) = -ln(rce/2(Dt)^*^-0.29 



Gg = Geothermal gradient 

Tge = Geothermal temperature at depth of fluid entry 

Tfe = Fluid entry temperature 

Z = distance from entry zone 

Q= Mass flow rate 

□f = Fluid density 

Cf = Fluid specific heat 

kh = formal thermal conductivity 

U = Overall heat transfer coefficient 

t = Flowing time 

Tci = Iimer radius of casing 

Tee = Outer radius of casing 

□ = Thermal difiusivity of casing 

Equivalent algorithms exist for calculating the flow rate of fluid being injected into an 
underground reservoir. See for example the paper '''Temperature Logging in Injection 
Wells '\ Witterholt E J and Tixier M P, SPE 4022. Similarly, for gas production, see for 
example the paper '^Temperature Surveys in Gas Producing Wells by Tixier M P and 
Kunz K S, AIME Annual Meeting, Chicago 1955. These algorithms may be improved 
by taking into account the effective heating due to flowing friction pressure drop. 
Additionally, changes in pressure, volume and temperature properties up the well may be 



10 



1 



taken into account by the use of computer nodal analysis. Suitable commercially 
available temperature modelling software includes that sold by Landmark under the trade 
name WellCat. 

5 During the calibration period, the algorithms are used while the actual fluid flow rate is 
determined independently, so that the flow rate becomes a known value. The known 
flow rate, whether a flow rate of fluid produced from the well or a flow rate of fluid 
injected into the well, is compared with the measured temperature profile of the length of 
conduit in consideration, as a function of producing time and depth. The comparison 

1 0 may be optimised by modification of the formula constants such as fluid specific heat 
(Cf), the thermal conductivity of the surrounding rock formation (kh), the overall heat 
transfer coefficient (U) and the thermal dififusivity of the casing (□), which are specific 
to a particular well completion, a particular formation, and particular fluid properties. In 
effect, the calibration phase in which the actual flow rates are knovra enables the 

1 5 constants in the algorithms to be established for a particular well. A more accurate 

analysis may be obtained by calculating the formula constants themselves, by employing 
a least squares regression fit of the predicted data to the measured data as a fiinction of 
time and depth and flow rate. 

20 Once the constants in the algorithms have been derived as a function of known flow rates 
during the calibration period, it is possible to calculate the flow rate in an observation 
period using the same derived constants. Normally, the observation period will follow 
the calibration period, and this is of course essential if real time data is required. It 
would however be possible to use an observation period before a calibration period, if it 

25 would be acceptable to derive only historic data. 
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A second preferred method for deriving mass flow rates from temperature data occurs 
where the temperature of the flowing fluid changes as a result of a change, typically a 
loss, of pressure. This may be due to a discontinuous change in the size or type of 
5 conduit, for example a 3/4" to 4/4 " (89mm to 1 14mm) crossover, or at a tubing shoe or 
the like, or a loss of pressure along a horizontal length of production conduit. The 
temperature changes can be related, by the use of appropriate equations, to the flow rate 
of the fluid. Again, parameters that are fixed as constants in a particular downhole 
environment are the thermal properties of the flovraig fluid, of the installed completion 
1 0 and of the surrounding formation. Dynamic properties of the fluid such as its Joule- 
Thompson and PVT characteristics should also be taken into account. 

When this approach is followed, the temperature profile at a change in flowing cross 
section of the conduit is characterised by either an increase or a decrease in temperature. 
15 The mass flow rate of the fluid can be determined, as before, from an analysis of the time 
based distributed temperature profile over the length of conduit concerned, where the 
PVT characteristics of the flowing fluid are known and the other relevant constants are 
derived from measurements made during the calibration period when the flow rates are 
known., . 

< 

20 

The algorithms defining the relationship between the temperature profile, the length of 
the conduit being investigated, and time, with respect to the thermal properties of the 
flowing fluid and the surroundings, are available from published literature describing 
heat transfer in pipelines. Reference is made in particular to two books: 

25 
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1 Hein, Michael A. "HP 41 Pipeline Hydraulics and Heat Transfer Programs 
PennweU books, 1984, ISDN 0878 14 255X 

2 Carslaw, H.S. and Jaeger, J.C. '''Conduction of Heat in Solids " Clarendon Press, 
5 2°^ Edition, 1959. 

The invention is illustrated diagrammatically by way of example in the accompanying 
drawing, Figure 1 , which illustrates on the left hand side a cross section of a producing 
oil well, and on the right hand side a graph of temperature (abscissa) against height/depth 
10 (ordinate), the depth scale being indicated by the corresponding location in the adjacent 
depicted well bore. 

In Figure 1 , the oil well is shown with a casing 1 1 extending from the ground surface 1 2 
into and through a producing reservoir 13. Production tubing 15 extends inside the 
15 casing from the usual oil flow control apparatus 16 located above ground at the wellhead 
and terminates inside the casing at a depth D above a producing zone 17. The upper 
boundary of the producing zone is marked by a closure 1 8 which holds the lower end of 
the production tubing in place within the casing. 

* 

20 In accordance with the invention, an optical fibre is deployed within the well in a 
suitable duct, such £is optical fibre deployment tube 20, which is a continuous tube 
having two limbs, in effect a U-tube, beginning and ending in connection with surface 
mounted instrumentation 22, including a light source, a detector and data processing 
apparatus, as means for passing light along the optical fibre, receiving hght emergent 

25 therefrom and interpreting temperature- and location-related characteristics of said 
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emergent light in terms of the temperature profile of the optical fibre at a series of 
locations along the fibre sensor deployed in the well. The deployment tube extends firom 
the instrumentation down the well between the production tubing and the casing, through 
closure 1 8 and through the producing zone 1 7, returning by the same route. The tube is 
thermally conductive, and will typically be clamped to the outside of the production 
tubing, in good thermal contact therewith, but may alternatively be installed on the oil 
well casing, if fluid temperatures in the annulus between tubing and casing are to be 
measured. 

The instrumentation 22 may comprise commercially available instrumentation, such as 
the model DTS-800 made by York Sensors Ltd of Chandlers Ford, Hampshire. The 
optical fibre within tube 20 is desirably coated to withstand the high temperatures and 
corrosive fluids encountered in a downhole environment. Typically, the fibre will 
comprise an inner core or wave guide of about 50|im diameter surrounded by a lower 
refiractive index cladding, total diameter about 125|im. ,The cladding may be coated with 
a sealing layer that is impervious to downhole fluids, and finally an abrasion resistant 
coating to bring the total diameter of the fibre up to about 155-400iim. 

Such a fibre can be deployed in the deployment tube 20 by hydraulic means, and can 
correspondingly be replaced if necessary. The deployment fluid provides a thermal 
bridge to the tube wall, so that the optical fibre is in thermal contact with the fluids 
whose mass flow rates are to be measured, in the conduit comprised of the production 
tubing and the well casing. With appropriate instrumentation, temperatures can be 
measured firom -40DC to +300DC, with an accuracy of 0.5D or better. The temperature 
measurements are absolute rather than relative, do not drift and do not need calibration 

14 



against any reference. In practical terms, the information can be obtained continuously, 
meaning as frequently as every 7 seconds at intervals of Im along the fibre, for a fibre 
length of up to 10km. In a 40km fibre, temperature readings at 10m intervals are 
effectively continuous over the length of the fibre. * 

The graphical right hand side of Figure 1 illustrates the measured temperature profiles. 
The natural geo thermal profile 30 is a straight line relation between depth and 
temperature. It is derived from simple temperature measurements, by conventional 
apparatus, at different depths. Deviations from a straight line, due to varying thermal 
conductivities in the surrounding formation, are considered to be so minor as to be 
insignificant for the purposes of the invention, in almost all cases. 

The heavier curve 32 represents the distributed temperature profile over the whole well. 
It coiQcides with the geothermal profile below the reservoir 13, where there is no flow, 
and the fluid is at equilibrium with its surroundings. As fluid enters the well from the 
reservoir and rises in the producing zone 17, it passes into cooler regions and begins to 
lose heat to the surrounding formations, which act as a heat sink at a temperature related 
to depth by the geothermal profile. Depending on flow rates, conductivities and the like, 
the temper-ature of the fluid rising in the well falls. 

m 

At depth D, the diameter of the conduit, constituted by the casing 1 1 in the producing 
zone 17 and by the tubing 15 above depth D, is suddenly reduced. As illustrated, the 
fluid temperature drops sharply at poiut 36 as it enters the narrower bore of the 
production tubiag. Thereafl:er, the rising fluid contiaues to show a temperature 
differential to the surrounding formation, and the relation between its temperature and its 
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depth approaches the asymptote 34 which is a straight line at the geotherraal gradient, 
displaced from the geothennal profile 30 by the difference in temperature between the 
temperature of fluid rising in a steady state from an indefinitely deep well and the 
temperature of the surrounding rock. 

The distributed temperature profile curve 32 provides the basis for the derivation of the 
required mass flow rates, as described. It is calibrated by measurements made at a time 
when mass flow rates are determined independently by conventional means, such as by 
spinner or venturi methods. Data processing in the instrumentation 22 provides real 
time information which is invaluable in managing the oil well. The novel application of 
thermal analytical techniques to temperature data derived from optical fibre distributed 
temperature sensors, in accordance with this invention, enables accurate, substantially 
non-intrusive, and easily replaceable downhole apparatus to give continuous real-time 
mass flow data. This in tum can be used in many ways, as well known in the art, to 
enhance oil well management. 
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